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ASTEROID-TYPE ORBIT EVOLUTION NEAR THE 5:2 RESONANCE." 

S.I. Ipatov, j jJ 

M.V. Keldysh Institute of Applied Mathematics, Moscow, USSR f 


In the case of the 5:2 commensurability with the motion of 
Jupiter an asteroid can reach the orbits of Mars, Earth and Venus 
when eccentricity e is greater than 0.41, 0.65 and 0.74, 

respectively. Tor individual fictitious asteroids Ipatov [1] and 
Yoshikawa [6J obtained a growth in e from 0.15 to 0.74-0.7o. Rates 

of changes in orbital orientations are different for Mars, Earth, 
Venus and the asteroid. Therefore, for corresponding values of e, 

the asteroid could encounter these planets and leave the gap at 
those encounters. In order to investigate this hypothesis of the 
5:2 Kirkwood gap formation, Ipatov [2] studied the regions of initi- 
al data for which the eccentricities of asteroids located near the 
5:2 commensurability exceeded 0.41 during evolution. The orbit evo- 
lution for 500 fictitious asteroids was investigated by numerical 

integration of the complete (unaveraged) equations of motion for 

the three-body problem (Sun-Jupiter-asteroid). The equations of 

motion were integrated in the time intervals .7^5 *10 J tj (tj is the 

heliocentric orbital period of Jupiter) in the planar model,,. 

7>,10 4 tj at initial inclination 5°^/ o ^20° and 7=10 5 tj at i* o =40°. The 

larger interval T was taken at * o =40° because in this case for the 

majority .of runs maximum values of e and i were reached in the time 
At>2-10 4 tj. 

Various initial orientations of orbit of the asteroid and its 
location in orbit were considered when initial value of asteroidal 
semimajor axis a Q was equal to the resonance value a 5/2 . It was 

obtained that maximum value e „„ of asteroidal eccentricity during 

evolution exceeded 0.41 and 0.65, respectively, for 2/3 and 1/3 of 

all investigated asteroids with the initial asteroidal eccentricity 

e o =0.15, / Q ^20 °, a cT a 5/2 an< ^ th e present-day value of Jupiter’s 

eccentricity ej . If at a Q =a 5/2 anc * some starting values of e, i, 

the argument of perihelion w, the longitude of the ascending node Q 
and the true anomaly v it was obtained that e max -0-41, then with 

the same starting values of e, i, w, £2 and v we have e max -0-41 also 

for the values a Q located in some vicinity of a g/2 . For all the 

fictitious asteroids examined with / =40° we obtained e rTlov ^0.6. 

Since earlier ej exceeded 0.06, we took ej= 0.06 when determining 

the maximum region of a Q and e Q for which fictitious asteroids with 

some initial orbital orientations were Mars crossers. It was shown 
that the outer boundaries of this region coincided wittL the 
boundaries of the 5:2 Kirkwood gap. For £ Q -0.2 and < o ^20° the 

regions of initial data for which fictitious asteroids were Earth 
crossers and Mars crossers are close to each other. Since the 
radius of Earth is approximately twice that of Mars, it is likely 
that meteorites, after having migrated from the 5:2 gap, impacted 
Earth more oiten than Mars. Both the range free of real asteroids 
for which £>5° and the range of a Q and e Q with 5°^' o ^20° for which 

e max~ 0 ' 41 are larger to the right of the resonance (for a>a 5/2 ) 
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than to the left (particularly for e^0.2). 

Analytical studies of the circular three-body (Sun-Jupiter- 

asteroid) problem, carried^ out, by Williams [5], ^howed that e 
reached a maximum when Aw=w— Wy=0 or Aw=180 u where w is the longitu- 
de of asteroidal perihelion and its value for Jupiter is designated 
by subscript This relation between e and Aw was obtained for 

most of the fictitious asteroids. Different cases of the interrela- 
tion of the time variations in e and Aw have been examined 
repeatedly [1-3, 6-7J. Our calculations showed that asteroids are 

Mars crossers and Earth crossers as a rule for certain types of 
interrelations of variations „in e and Aw. For the most of consider- 
ed fictitious asteroids the osculating elements changed in time 

almost periodically. We can distinguish several periodical compo- 
nents in the plotted time dependence of e and Aw with different 
amplitudes which may vary during evolution. Let us denote by T £ the 

period of variation in e with the largest amplitude, and by F~ the 

period of long-term variations in Aw if Aw librate around some 
constant usually equaj^ to 0 or 180°, or the time during which Aw 
change by 360® if Aw circulate. It was obtained for most runs that 

7~=7\,«const. In this case we can distinguish several types of 

interrelations N~ for variations in Aw and e. They are characteri- 
zed by the formula Aw=360 o (f-/ + )S~/T e +£+Aw£ and by the values of 
S~ , Awt , <S~ , Aw~ , S~ and N n presented in Table T. We have N =R 
if a value of asteroidal semimajor axis reaches a 5/2 during 

evolution and N a =N if it doesn’t reach a 5/2 - The first term in the 

right-hand ^side of this formula is the monotonic component in vari- 
ation of Aw, with Sgj denoting the sign of this component. There is 

no monotonic component when S~=0. The term Aw* denotes the value of 

Aw at time t=t + when e (t+)~ e max * an< * Aw~ is the value of Aw at 

time t=t_ when e(0~ e m / n • The £ is the oscillatory component in 

the variation of Aw. It was calculated that during evolution 
|£|<180° and £-0 for t=t + ±k*T £ /2 , where k=0 , 1, 2, 3, ... We 

introduce the plot of w(/) with w(f) averaged over the shot-period- 

ic oscillations. It was considered that S~=l if 3w/3t>0 at t=t , and 

U> T 

S~=-l if 3w/dt<0. The values of S~ for t=t are denoted in the same 
w w — 

way. - For the value of Aw almost didn’t changed with time and 

it was considered that S~=S~=0. When Aw=0 for N~=B or N~=H the 

w w w w 

eccentricity has, as a rule, a local maximum which is several times 

lower t han the main maximum obtained at Aw=180 . This local maximum 

is higher than the minimum by 0.05-0.1 for and by 0.02-0.05 

for N~=H. 

In some cases at / o =40° the eccentricity doesn’t reach maximum 
at Aw=0 or Aw=180°. In these cases some types shown in Table TT 

strongly depend on N a , Sq and S^. If Aw=w-Wjr librated during 
evolution about 0 or 180° with the amplitude du&20°, then it was 
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assumed that 5 =0. For all other investigated cases the Aw 

W 

increases during evolution (neglecting short-periodical variations) 
and 5=1. The variable R decreased during evolution if k 90 ana 

increased if £>90°. The values of show the sign of changes in 

R. For all types N~ presented in Table T it was obtained that 5^=1 

and 5 q=- 1. The type jV~=y t is a particular case of the type N~=y. 

Subtypes Af~=t/ - and N~=U~ h , N~=W~ and N~=W + are the parts of types 

N~=U and N~=tF, respectively. Each of these subtypes is replaced by 

U/ W a 

other subtype when i reaches 90 . It was obtained for types N~=U 
and that e max > 0-99 and the maximum value of i during evolution 

<max~160°' We obtained for one asteroid that e max = 0 .9999. Types /, 
T, E, Y t as well as all the types presented in Table TT were not 
obtained by Yoshikawa [7]. 

TABLE I 

Types of interrelations N~ for the long-term variations in e and Aw 


N » 

B 

c 

D 

H 

J 

Q 

T Z 

0 

I 

E 

y , 

N a 

R 

R 

R 

R 

R 

R 

R R 

N 

N 

N 

N 

5 w 

1 

0 

-1 

0 

0 

1 

0 -1 

1 

0 

0 

-1 

Aw+ 

180* 

3 180 c 

* 0 

0 

0 

0 

O 

o 

00 

rH 

o 

o 

00 

V-* 

0 

0 

90° or 180° 

180 ( 

4 

1 

1 

-1 

-1 

1 

1 

0 -1 

1 

1 

1 

-1 

Aw - 

e 

0 

Oo 

o 

o 

00 

o 

o 

o 

0 

00 

o 

00 

o 

o 

o 

h-* 

00 

o 

o 

0 

90° or 180° 

0 

s~ 

1 

-1 

-1 

1 

-1 

1 

0 -1 

1 

-1 

-1 

-1 


TABLE n 

Types N~ when e does not reach a maximum at Aw=0 or Aw=180°. 


N~ 

w 

X 

w~ 

W + 

Y 

u~ 

U + 

N a 

R 

R 

R 

N 

N 

N 

s a 

-1 

-1 

1 

-1 

-1 

1 


0 

1 

1 

0 

1 

1 


Let us consider interrelations of variations in orbital 
elements with periods less than T for the three-dimensional model. 

The times during which Aw and R change by 360° are denoted by T 

and Tq . For the types N~=X and N~=Y the value of T ' is the period 

of libration of Aw. Periods of long-term variations in i are 

denoted by T- and i. For most of the considered asteroids 

T^=T . Computer simulation results showed that for variations in i 

(r C 

with the period T- the local maximum of e always corresponds to the 

minimum of i and the local minimum of e - to the maximum of i. It 
was obtained that 27=7^ for most computer runs with / Q *0. For this 
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equality the minimum of i was reached at Aw=0 or at Aa)=180 o . The 

greater are the values of e Q and i Q , the greater are the 

variations in i and e with period T^T^/2. Two other equalities 

between T. and T ' : T-=T f . and T—2T.. were obtained more rarely. For 

l w l w l iv 

both cases i reached minimum at Aa>=180 . 

The relation 7 ^=Tq was fulfilled as a rule; besides for the 

considered runs i reached minimum at £2=180° and maximum at Q=0 for 
those cases when we observed no variations in e and i with the 

period T—T./2. Interrelations of periods of variations for four 

orbital elements (/, e, Q and o>) were obtained mainly for 5 -* o -10 

and e Q ^0.15. For some asteroids with Af~=/ the relation 

T*.=T~T,=T =Trt=T~ was fulfilled and in this case the minimum i and 

the maximum e were reached at 0=180 and Aa>=180 , while the maximum 
i and the minimum e — at £2=0 and Aco=0. If Af~=0 , it was obtained 

for most of the considered asteroids that T --T =2F \j=T e =T ^ when 

*0=10° and e Q ^0.05. 

For some fictitious asteroids the limits A a of variation in a 
exceeded the width of the gap and the asteroids could migrate from 

one side of the gap to the other side. Maximum values of A a/aj for 

types N~ denoted by C, T, W and U reached 0.01 and for N~=B and 

N~=Q - 0.013 and 0.016 respectively. 

In the three-dimensional space ( e max > a Q * e G )» the surface 

of values of e ov has a number of "plateaus" separated by steep 
m3 x 

"cliffs". Each plateau corresponds to a certain type N~ . The 

cliffs were obtained at the boundaries of initial data 
corresponding to different types N~ (excluding the transitions 

between types N~=Q and N~=0 as well as A/~=/ and N~=I) and were 

almost absent in a number of cases with e o =0.3. With initial data 

near these boundaries the transitions between the types were 
obtained for six asteroids in the planar model and for twenty 

asteroids in the three-dimensional model. Sidlichovsky [4] obtained 
that an orbit was chaotic if initially small starting eccentricity 

exceeded 0.3 during evolution. At types N~ denoted by B, C, D, H, 


W, Z, E and V for all considered asteroids it was obtained that 

e >0.41 and be=e v -e mi „-0.24. Therefore, the number of chaotic 

max max min 

orbits may be considerably larger than the number of asteroids for 
which such transitions were obtained. 
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